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Abstract 



We consider symmetric binary mixtures consisting of spherical particles 
with equal diameters interacting via a hard-core plus attractive tail potential 
with strengths e^-, i,j = 1,2, such that en = £22 > £12- The phase diagram 
of the system at all densities and concentrations is investigated as a function 
of the unlike-to-like interaction ratio 5 = £12 /en by means of the hierarchical 
reference theory (HRT). The results are related to those of previous investiga- 
tions performed at equimolar concentration, as well as to the topology of the 
mean- field critical lines. As 5 is increased in the interval < 5 < 1, we find 
first a regime where the phase diagram at equal species concentration displays 
a tricritical point, then one where both a tricritical and a liquid-vapor critical 
point are present. We did not find any clear evidence of the critical endpoint 
topology predicted by mean-field theory as 5 approaches 1, at least up to 
5 = 0.8, which is the largest value of 5 investigated here. Particular atten- 
tion was paid to the description of the critical-plus-tricritical point regime in 
the whole density-concentration plane. In this situation, the phase diagram 
shows, in a certain temperature interval, a coexistence region that encloses 
an island of homogeneous, one-phase fluid. 
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I. INTRODUCTION 



A major difference between the phase behavior of one-component fluids and binary 
mixtures is that, even if one considers just simple systems with a spherically symmetric, 
Lennard- Jones (LJ) like interaction profile, the qualitative features of the phase diagram 
of mixtures depend very sensitively on the parameters of the microscopic potential. If the 
interaction between a particle of species i and a particle of species j is modeled as the sum 
of a hard-core repulsion and a longer-ranged attractive tail, the relevant parameters are the 
hard-sphere diameters and the strengths of the attractive contributions, i,j = 1,2. 
This parameter space is drastically reduced by focusing on a particular class of systems, 
generally referred to as symmetric mixtures, such that an = 022 = <Jvi — e n = e 22 = e- 
Since the quantities o and e can be included into the definition of the temperature T and 
number densities pi by introducing standard reduced units p* = Pier 3 , T* = kftT/e, being 
the Boltzmann constant, it follows that the only parameter affecting the phase diagram is 
the ratio of the interaction strengths between unlike and like species 5 = e 12 /e. 

Clearly, symmetric mixtures appear quite artificial when considered as a model of real 
binary fluids. In fact, some of the features of their phase behavior hinge on the invariance 
with respect to the exchange of the two species, and are not found in the phase diagram 
of real mixtures. As has already been observed [|T| , symmetric mixtures are better seen as 
a model for a one-component fluid, whose particles have been endowed with a two-state, 
spin-like variable in addition to their translational degrees of freedom, so that their mutual 
interaction depends both on their relative position and on their "internal" state, namely 
whether the interacting particles belong to the same species or not. As such, this model 
mixture is closely related to other models of dipolar ||[| and magnetic fluids, 
especially Ising-spin fluids [^pyiQf . In these systems, the phase behavior results from the 
interplay between the liquid-vapor phase separation and the additional transition, e.g. para- 
ferromagnetic, associated to the spin-like degrees of freedom. In symmetric mixtures, the 
latter corresponds to the mixing-demixing transition. 

Because of the relative simplicity of this model compared to a generic binary mixture and 
of the possibility of generating the whole spectrum of phase diagrams by acting on just one 
parameter, symmetric mixtures have been widely studied both by mean-field theory [|l| and 
by numerical simulations []l^ [lTr| , [l^ , |T^| . p~^| . |T3| , p~0| | . A situation which has been given special 
attention is that of equal species concentration x = p2li.P1 + P2) = 1/2. In this case, 
accurate numerical simulations fl|] have qualitatively confirmed the mean-field scenario for 
the phase diagram as the parameter 5 is varied in the interval < 5 < 1. This will be 
considered in detail in Sees. |T|, [IV|. Here we just recall that the mean-field phase diagram at 
equal species concentration presents both a first-order coexistence boundary which separates 
a low-density fluid from a high-density one, and a line (the so-called A-line) of mixing- 
demixing critical points. Beyond the A-line, the fluid actually consists of two demixed fluids 
in equal amounts, one at a certain concentration x, and the other at concentration 1 — x, 
so that the overall concentration of the two species remains the same. For large enough 5 
(Si < S < 1 with 5i = 0.708 according to the mean-field result [0) the coexistence curve 
ends into a liquid- vapor critical point, while the A-line intersects the coexistence curve at 
a point of first-order transition, thereby terminating into a critical endpoint. At small 5 
(0 < 5 < 5 2 with 62 = 0.605 in mean field |l|), on the other hand, the point at which 
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the coexistence curve meets the A-line coincides with its critical point. The latter is then 
referred to as a tricritical point, since on approaching this point from low temperatures, one 
observes the simultaneous coalescence of three phases, namely the low-density vapor and 
the two demixed high-density fluids. Finally, in a narrow interval of S values 62 < S < Si 
intermediate between those corresponding to the two topologies described above, one has the 
occurrence of both a liquid-vapor critical point and a tricritical point. Despite the qualitative 
agreement between the mean-field scenario and the simulation results, there are still several 
points that deserve further investigation. First, mean- field theory and simulations show 
considerable quantitative discrepancies, which concern both the position of the critical loci 
and the values of S at which the changes in the topology of the phase diagram occur. This 
is in itself not surprising, as mean-field theory cannot be expected to be quantitatively very 
accurate. Therefore, one would like to go beyond it by means of a theoretical treatment 
which includes fluctuations in the order parameter of the transition, be it of the liquid- vapor 
or mixing-demixing kind. Two relevant issues in this respect are whether the mean-field 
scenario is qualitatively recovered even after fluctuations have been taken into account, and 
which is the extent of the quantitative changes involved. Moreover, the case of equal species 
concentration corresponds just to a certain plane, albeit undoubtedly of special interest, of 
the space of thermodynamic states. Simulation studies of the phase diagram have indeed 



been performed also at fixed density p and variable concentration x p4yi6| , but mapping 
the phase diagram in the whole thermodynamic space for different values of the parameter 
S would require an exhorbitant number of simulation runs, and therefore hardly appears as 
a viable strategy in view of the computer time required. It is then tempting to resort to 
theory in order to explore the phase diagram at general density and concentration and find 
out how it changes by changing S, so as to see what the phase diagrams corresponding to 
the three regimes outlined above look like, as one moves away from the x = 1/2 plane. 

We aim to address these topics in the present work by means of the hierarchical reference 
theory (HRT) of binary fluids. This theory has already been successfully applied to the 



description of the critical behavior [ I7|| , crossover phenomena (13113 > anc ^ P nase diagram 1 19 



of simple fluid mixtures. Few results for the phase behavior of symmetric mixtures have 
also been reported fll9| , |20|| , but a systematic investigation of this system by HRT has not 



been undertaken yet. In our opinion, HRT is especially well suited for such a study, with 
particular regard for the tasks stated above. In fact, the aim of this approach is determining 
how the Helmholtz free energy of the mixture is affected by the introduction of density 
and concentration fluctuations. This is achieved via a renormalization-group (RG) like 
procedure, where the long-wavelength Fourier components of the microscopic interaction 
are gradually introduced into the hamiltonian of the mixture. Any intermediate stage of 
this process, such that only Fourier components with wavevectors exceeding a certain cutoff 
Q have been taken into account, physically corresponds to suppressing fluctuations on a 
lengthscale L > l/Q. Long-range fluctuations are recovered in the limit Q — > 0, when 
the free energy of the fully interacting system is obtained, while the mean-field free energy 
enters as the initial condition at Q = 00. The main advantage of HRT over other liquid-state 
theories is that it embodies several features of the RG description of critical phenomena in 
a treatment based on the microscopic hamiltonian of the fluid. These include scaling, non- 
trivial critical exponents, and the correspondence between universality classes and different 
fixed points of the RG flow [TTj . A fact of particular relevance for the investigation pursued 
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here is that the inclusion of long-range fluctuations has the effect of preserving the correct 
convexity of the free energy in the whole thermodynamic space. Whenever phase coexistence 
occurs, one does not find any domain of instability as in the mean-field approximation, and 
the conditions of thermodynamic equilibrium between the phases at coexistence are enforced 
by the theory itself. At each given temperature, the coexistence region is then immediately 
recovered as the locus in the density-concentration plane where the chemical potential of 
each component is constant along the lines of fixed pressure, with no need of imposing this 
condition a posteriori by a Maxwell construction. For binary mixtures, the latter proves to be 
quite cumbersome already at the mean-field level, and is much more so for more sophisticated 
integral-equation theories, in which the occurrence of phase separation generally entails the 
presence of some forbidden domain, where the theory cannot be solved at all. Therefore, the 
ability of straightforwardly mapping the phase diagram is a valuable asset of HRT. This is 
especially true in the present case where the topology of the phase diagram is very sensitive 
to changes in S, while at the same time, as will be seen in the following, the features that 
allow one to discriminate between different topologies are often detectable only in a narrow 
window of the thermodynamic space. Because of the lack of a solution defined for every 
(T, p, x) state, pinpointing all these features by conventional integral-equation theories would 
undoubtedly prove extremely difficult, perhaps even impossible. 

In this work we have considered symmetric mixtures of additive hard spheres interacting 
via an attractive Yukawa tail potential Wij(r) = — crej ? -e~^ r / fT ~ 1 )/'>'', where r is the interparticle 
distance and z is the inverse range of the interaction, which has been fixed to the value z = 1.8 
for both like and unlike species. This hard-core Yukawa (HCY) form has been preferred to 
the square- well potential used in Ref. M. The latter lends itself well to simulation, but the 
very slowly decaying behavior of its Fourier transform makes it somewhat tedious to use in 
HRT. The HCY potential has already been adopted in a number of studies of symmetric 
mixtures based on the mean spherical approximation (MSA) plj , the optimized random 
phase approximation (ORPA) ||22fl , and the self-consistent Ornstein-Zernike approximation 
(SCOZA) p3| , all of which yielded for the phase diagram at equimolar concentration the 
same behavior found in mean-field theory. The scenario that comes out of our investigation 
by HRT agrees qualitatively with the mean- field one in predicting that, as 5 is increased, the 
phase diagram at equimolar concentration exhibits first a tricritical point, and subsequentely 
both a tricritical and a liquid- vapor critical point. However, we did not find any clear 
evidence of the occurrence of the mixing- demixing critical endpoint given by mean-field 
theory as 5 approaches 1, at least up to 5 = 0.8, above which further investigation is hindered 
by the finite resolution of the density grid used in our numerical calculation. Besides the 
scenario sketched above at equal species concentration, other interesting features of the 
phase diagram emerge as we move off the x — 1/2 axis. For instance, the intermediate 
regime where both a critical and a tricritical point are found for x — 1/2 is marked in the 
p-x plane by a coexistence region which in a certain temperature interval contains a "hole" 
or "island" of homogeneous, mixed fluid. Double critical points and, for a certain value of 
S, tricritical points are also observed for unequal species concentration. 

The paper is structured as follows: the HRT for a binary fluid is described in Sec. [H|. 
Mean-field theory is recovered within this approach as a zeroth-order approximation, and 
the mean-field results for the critical lines of a symmetric mixture as the parameter 5 is 
varied in the interval < 5 < 1 are shown in Sec. |TT|. The HRT phase diagram for different 
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values of S is discussed and compared with the mean-field predictions in Sec. [IV]. Finally, in 
Sec. we summarize our findings and draw our conclusions. 



II. HRT EQUATIONS 

Here we briefly review the HRT approach for a binary fluid. A more detailed derivation 
can be found in previous works |T7| , |I9| , pi| . 

We consider a model mixture consisting of particles of two species interacting via a two- 
body spherically symmetric potential Vij(r), where the indices i, j label the particle species. 
The derivations presented in this Section do not hinge on the fact that one is dealing with 
a symmetric system such that t>n(r) = f22( r ) ; so they are equally valid for non-symmetric 
systems. We assume that %(r) can be split as the sum of a singular contribution fJHr) 
which accounts for the short-range repulsion between the particles, and a longer-ranged, 
smooth attractive tail Wij(r) which may induce fluid-fluid phase separation. The properties 
of the mixture interacting via the repulsive potential vfj(r) alone are considered as known, 
so that it acts as a "reference" or unperturbed system. For a hard-core plus tail potential, 
the natural choice for the reference system is just a binary mixture of additive hard spheres, 
which can be described by the Mansoori-Carnahan-Starling-Leland equation of state 



and the corresponding Verlet-Weis parameterization for the two-body correlations p6|,|2T 



We remark that here we will not be concerned with a possible demixing transition occurring 
in the hard-sphere mixture, since such a transition may come along only as a consequence 
of depletion interactions when the particles differ widely in size. In the present case of eq- 
uisized particles, the reference system reduces to a one-component hard-sphere fluid, and 
all the fluid-fluid transitions displayed by the system are necessarily due to the attractive 
perturbation Wij(r). The HRT differs from the conventional liquid-state approaches in the 
way this perturbation is dealt with. In order to accurately describe the long-range fluc- 
tuations that are important in criticality and phase separation, the attractive part of the 
interaction is switched on gradually by introducing a Q- system with a modified interaction 
v ij( r ) = v fj( r ) + w ?j( r )> where w®{r) is defined in such a way that its Fourier transform 
wfj(k) coincides with that of the original attractive potential Wij(k) for k > Q, and vanishes 
for k < Q. Inspection of the diagrammatic series of the Helmholtz free energy of the mix- 
ture in terms of the perturbation Wij(r) and the correlation functions of the reference system 
shows that introducing such an infra-red cutoff in the interaction is physically equivalent 
to inhibiting fluctuations with characteristic lengths L > 1/Q. If Q is made evolve from 
Q = oo, the Q-systems evolve from the reference system by acquiring fluctuations of longer 
and longer wavelengths. The fully interacting system is recovered as the Q — > limit of 
such a process. Only in this limit true long-range correlations are allowed to develop in 
the fluid. The equation for the corresponding evolution of the Helmholtz free energy Aq 
of the Q-systems can be determined exactly and is related to the attractive perturbation 
in momentum space &ij(k) = —/3wij(k) where (3 = l/(k^T), and to the direct correlation 
function of the Q-system in momentum space cfj(k). We recall that the direct correlation 
function is related to the structure factor of the fluid by the Ornstein-Zernike relation. If 
we denote by cq(/c) the 2x2 symmetric matrix with elements cfj(k), for a binary mixture 
this relation reads 
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(1) 



where Cq{k) is the inverse matrix of cq(k) and SfAk) is the partial structure factor of the 
Q-system. Here, unlike in the convention commonly adopted in liquid-state theory, cfj(k) 
contains its ideal-gas contribution —5ij/pi. The evolution equation for the Helmholtz free 
energy is most easily formulated in terms of a modified free energy j^q and direct correlation 
function 'tfn(k) defined as: 



«#(*) = $(A0 + -*§(*). 
The HRT equation for s/q is then 



1 2 



d 3 k 

(2tt) 3 



<9g 



Q 2 



4tt 2 



io g {det [i-^ 9 1 (g)$(g)]} 



*«(*)- SjftAOJ (2) 
(3) 



(4) 



where again <&(Q) and I^q (Q) are 2x2 symmetric matrices, the latter being the inverse of 
the matrix ^q(Q). We note that for Q — ► 0, i.e. at the end of the evolution process, 
and $(&) coincide, so that in this limit the modified quantities sfq, yield respectively the 
true free energy and direct correlation function of the fully interacting system. For Q = oo 
instead one has $q(A;) = 0, and s#q and are nothing but the mean-field free energy 
and the random-phase approximation (RPA) direct correlation function in the presence 
of the full perturbing potential $(&). These play the role of the initial conditions of the 
evolution equation (^), which then describes how the mean-field estimate for the free energy 
is affected by the inclusion of fluctuations. This equation is manifestly not closed, since the 
evolution of the free energy $$q is related to the matrix of the direct correlation function 
c tfQ(k), which is itself unknown. In fact, Eq. (|4|) is just the first equation of an infinite 
hierarchy for the direct correlation functions of increasing order: for instance, the evolution 



of is related to the 3- and 4-body direct correlation functions in Fourier space p4 



A point of crucial importance in order to implement a viable HRT scheme consists then 
in supplementing Eq. (Q) with some closure relation involving c €q{]z). Here, as well as in 
the previous applications of HRT, we have not resorted to the higher-order equations of the 
hierarchy. Instead, we have adopted for ^q(/c) an approximate form inspired by standard 
perturbative liquid-state theories: 



-■ij 



(k) + 1% $v(k) 



(5) 



where c^ s (k) is the Fourier transform of the partial direct correlation function of the hard- 
sphere reference system, which has been represented by the above-mentioned Verlet-Weis 
parameterization. The functional form of Eq. (|5|) for the direct correlation function is similar 
to that of the widely used RPA, which is recovered for i/Q — 1. In particular, both of them 
rely on the Ornstein-Zernike ansatz, i.e., the direct correlation function has always the same 
range as the potential, so that ^(k) is always analytic in k, including at the critical points 
of the system, where the real direct correlation function is instead expected to be non- 
analytic for Q — > 0. However, unlike in the RPA, the amplitude of the perturbation 
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is regarded as an unknown quantity, to be determined in such a way that each Q-system 
satisfies the compressibility sum rule. In a one-component fluid, this relates the isothermal 
compressibility to the zero-wavevector value of the structure factor, or equivalently of the 
direct correlation function in momentum space. Such a rule is readily generalized to binary 
systems, and in terms of the modified quantities #/q, ^(k) it reads 



By determining so that Eq. (|) is satisfied, one immediately finds from Eq. (|) 

<3Vn 



«£(*)= eg 8 (*) + 



dp l dp j 3 



<Pij(k) , (7) 



where we have set (pij(k) = $y(fc)/$y(A; = 0). Eq. ([?]) ensures that the Helmholtz free 
energy obtained by integration of Eq. @ is consistent with the compressibility route to the 
thermodynamics given by Eq. (|6]). This thermodynamic consistency condition plays a key 
role in the implementation of HRT. In fact, by substituting ^(k) as given by Eq. (0) into 
the HRT equation (|4]), a closed partial differential equation (PDE) for j^q is obtained, which 
involves both the first partial derivative of s^q with respect to Q, and its second partial 
derivatives with respect to the densities pi, p^. In order to integrate this equation numer- 
ically, we found it most convenient to cast it into a form where the partial derivatives of 
the unknown function appear only outside some "coefficients" that may depend both on the 
independent variables p\, p2, Q, and on the unknown function itself, but do not contain its 
derivatives. This allows us to take advantage of finite-difference schemes especially devised 
for equations of such a quasi-linear form, which combine robustness with a moderate com- 



putational cost p8fl . Both of these are very important requirements in our case. In fact, in 
order to deal with the divergence of the compressibility at criticality and phase coexistence, 
one has to resort to a very stable algorithm, while on the other hand for a diffusive PDE 
in three independent variables like the one considered here, the dimension of the vectors 
generated by the discretization procedure becomes rapidly very large even for a relatively 
coarse density step Ap. Hence, we have to use a solution scheme which is not too demanding 
computationally to prevent computer time from increasing beyond control. 

The method we adopted to re-write Eq. (§) supplemented by the closure relation (0) in 
quasi-linear form has already been illustrated in Ref. ||19|| , to which we refer the reader for 
details. Here we recall that, instead of s^q, we use as unknown function the quantity 

f/ = log{det[l-^ Q 1 (g)$(Q)]} . (8) 
This is found to satisfy a quasi-linear PDE of the form 

v dU T ,d 2 U T d 2 U nr d 2 U nT 

* U xi = K TA +L Wi> +M M + N (9) 

whose "coefficients" K, L, M, N, which will not be reported here, depend on the variables 
Pi; Pi-, Q both explicitly and implicitly via a set of three auxiliary variables. These can be 
identified either with the eigenvalues Ai, A2 of the symmetric matrix ^q(Q) and the angle 
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a of the rotation that casts &q(Q) into diagonal form, or with the corresponding quantities 
Ai, A 2 , 9 for the symmetric matrix of elements d 2 'j^q/ 'dpidpj. As discussed in [19|, the 



use of one or the other set of auxilary variables is dictated by the different behavior of the 
interaction in the high- and low-Q region. At high Q, where the Fourier transform of the 
attractive perturbation ifij(Q) typically displays oscillations, thereby vanishing for certain 
Q values, it is better to adopt as auxiliary quantities Ai, A2, 9. At small Q instead, when 
most of the attractive interaction has been included into the system and phase coexistence 
may occur, the set Ai, A 2 , a is the more appropriate choice. In both cases, the derivative of 
each of the auxilary variables with respect to Q can be expressed in terms of the variables 
themselves and the partial derivatives of U with respect to pi, p 2 . At any given Q, the 
resulting equations are used to update Ai, A 2 , a (or Ai, A 2 , 9), while the PDE (|) is used to 
update U. Since the matrix ^q{Q) that appears in Eq. (|8|) can be expressed either by Ai, 
A 2 , a or, via Eq. (|7|), Ai, A 2 , 9, it follows that U and the set of the three auxiliary variables 
are not mutually independent. The relation between U and the auxiliary variables is used 
throughout the integration procedure as a check of the accuracy of the calculation. 

In order to perform the numerical integration, we found it convenient to replace the 
independent variables pi, p 2 with the related variables £ = (pi + P2)c 3 , x = p\j (pi + p 2 ). At 
any given Q, U is then defined in the rectangular domain < 2 < 1, < £ < £ max , where 
the high-density boundary £ max was chosen equal to 1. At the beginning of the evolution 
process, i.e. for suitably large Q, the matrix ^q(/c) can be identified with the RPA expression 
for the direct correlation function which, as said above, is obtained from Eq. ([3D for Q = 00. 
This gives the initial condition for U via Eq. (H). The boundary conditions for U needed 
for the integration of Eq. (||) are determined as follows: for £ = the diagonal elements 
of the matrix ^q(Q) diverge because of the ideal-gas contribution to the direct correlation 
function — 5y/p» , so that U(£ — 0,x) vanishes identically. For x = and x = 1 only one 
of the two species is present. The corresponding boundary conditions are then given by 
the solution of the HRT equation for a pure fluid, which can be integrated numerically by 
specializing the procedure sketched above to a one-component system. The high-density 
boundary condition for £ = £ max is non-trivial because, unlike in pure fluids, we must allow 
the occurrence of phase transitions even at high density. As a consequence, we have to rule 
out the possibility of using for U at £ = £ max a simple approximation such as the RPA, as 
already done for the one-component case, because such a form would behave unphysically 
in the coexistence region. In general, we expect that at high density the compressibility 
of the fluid will be small, and concentration fluctuations will become more important than 
density ones. This corresponds to the transition becoming less liquid- vapor and more mixing- 
demixing in character. For the symmetric mixtures considered here, such an expectation 
holds rigorously: in fact, because of the special symmetry of the model, the high-density 
transition are of pure mixing-demixing type. As a consequence, we expect U to be much 
more sensitive to a change in concentration than in density. Accordingly, for £ = £ max the 
partial derivatives of U along £ in Eq. (|9|) have been disregarded with respect to those along 
the conjugate direction z — (pi — P2)c 3 - By switching from the variables £, z back to the 
variables £, x, we obtain a PDE in the variables Q, x on the boundary £ = £ max , whose 
solution yields the high-density boundary condition of Eq. (0). 

The numerical integration of Eq. (|9|) supplemented by the equations for the evolution 
of the auxilary variables was performed on a grid in the (£, x) plane containing 150 x 150 
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mesh points. The integration with respect to the variable Q was carried out by setting 
Q = Qoe~ f , t > 0, where the initial value of Q corresponding to t = was typically fixed 
at Qo = 30<7 _1 . The variable t was discretized using a step At = and the iteration in 
t went on until convergence in the quantity U was achieved outside the coexistence region. 
At low temperature this requirement can be satisfied for Q < 10 a . 



III. MEAN-FIELD CRITICAL LINES 

The phase diagrams of binary mixtures are usually classified according to the topology 
of their critical lines p9| , |30| . We then begin the discussion of the phase behavior of the 
symmetric mixtures as a function of the interaction parameter 5 by presenting the different 
shapes of the critical lines that are predicted by the mean-field approximation. We will then 
consider the HRT results for the phase diagram and elucidate the relationship between the 
critical loci and the behavior of the coexistence domains on changing the temperature. In 
doing so, we will also be in a position to compare the mean-field and HRT predictions. 

As observed above, the mean-field (MF) approximation enters the HRT as the initial 
condition at Q = oo, when no fluctuations have been introduced into the system: 

(3A MF pA m , 1 ^ A n ^ 1^ fd 3 k 



1 l r d k 

y -V=oc y + ^ E ftft%(^ ) - 2 Eft ' ( 10 ) 

The equations for the MF critical loci are also obtained within HRT as the lowest-order 
approximation to the requirement that the RG flow generated by the evolution equation (f|) 
drives the free energy of the mixture towards its fixed point. This gives the following 
equations [O: 




(11) 
(12) 

0, (13) 

where ipi, ip2 are obtained from the original densities pi, p2 via an a priori unknown rotation 
such that Eq. flT3| ) is satisfied. Eqs. (|TT|), ( [12] ) are formally similar to the equations for the 
critical point of a pure fluid, except that here the density p has been replaced by the variable 
ipi. We note that Eqs. flTT|), (|T3|) amount to requiring that the hessian determinant of the free 
energy must vanish at the critical point, and that the vanishing eigenvalue corresponds to 
the eigenvector directed along ipi. This is the linear combination of the densities that gives 
the direction of strongest fluctuation, and identifies the order parameter of the transition. 
For the solutions of Eqs. (|TTl) -(|l3|) to yield actual critical points, a further constraint has to 
be imposed, i.e. the condition of thermodynamic stability ensuring that even a critical point 
is an equilibrium state of the system, and as such is stable against density and concentration 
fluctuations. The stability conditions at a critical point read [|T7| 



9 




(14) 

1 d 2 f(3A MF \ «9 4 //M MF X 



<0. (15) 



The three equations (pTT|) (|X3|) with the conditions (0), ([15]) contain four unknowns, namely 
Pi, p 2 , T, and the state-dependent angle 7 of the rotation that identifies the proper axes 
■01, ip2- As 7 is varied, they will generate a set of lines in the thermodynamic space. We 
must observe that, when several critical lines are present, one should also check that a 
critical point does not fall into the coexistence region originating from another critical line. 
In mean-field theory, this circumstance can occur without violating the stability conditions 
(0)) (HI)' when the point considered lies between the binodal and the spinodal surfaces of 
a neighboring transition. In such a situation, a solution of Eqs. ( |11]) — (P~5|) corresponds to a 
critical point which, while still locally stable, is however globally metastable with respect to 
first-order phase separation. In order to assess this possibility, one should then determine 
the mean-field binodal surfaces by a Maxwell construction. This has not been done here. 
However, we can discriminate between stable and metastable regimes, at least for equimolar 
concentrations, by comparing our results with those obtained in Ref. |l[, where mean- field 
binodals were determined. We recall that at the mean-field level the phase behavior as a 
function of 5 is independent of the profile of the attractive interaction since this enters 
in the approximation only via its spatial integral. A change in the form of $jj is then taken 
into account by simply rescaling the temperature. 

Let us now consider the critical lines of the symmetric mixtures we are interested in. 
Because of the attractive interaction between the particles, we expect that for suitable 
temperature, density, and concentration, the system will undergo a liquid-vapor transition. 
This is certainly true for states at low or high concentration, where one of the two species will 
play a minor role. On the other hand, since the interaction between unlike particles is weaker 
than that between like particles, the internal energy will tend to promote demixing between 
the two species. For high enough density, this increase in the absolute value of the internal 
energy may overcome the loss in entropy implied by the demixing, and a mixing-demixing 
transition may appear. Therefore, both liquid-vapor and mixing-demixing critical lines are 
expected. As the unlike-to-like interaction ratio 5 is lowered, the mixing-demixing transition 
becomes more favored, and the corresponding critical line will move to lower density. The 
projections of the mean-field critical lines on the density-concentration plane for a relatively 
small value of 5 = 0.4 are reported in Fig. |I| The open dots mark the positions of the 
minima of the critical temperature, while the arrows drawn along the critical lines give 
a graphic representation of the relative weight of density and concentration fluctuations, 
thereby showing the direction of the order parameter ipi determined by Eqs. (p4])-(fT3|). 
Specifically, the angle <fi between the arrows and the density axis gives the fluctuation of the 
order parameter corresponding to a given fluctuation of the total density p and concentration 
x as 5ipi = Sp cos(f)+ p Sx smcf). Arrows parallel to the p- and x-axis then indicate pure liquid- 
vapor and mixing-demixing transitions, respectively. Fig. [I] shows a mixing-demixing critical 
line at equal species concentration x = 1/2. On the low-density side, this intersects another 
critical line that connects the critical points of the pure species. As one moves from the pure 
species to the equimolar mixture, the direction of the order parameter changes continuously 
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from pure liquid-vapor to pure mixing-demixing. The critical temperature, not shown in the 
figure (see the lower panel of Fig. |9] for the similar case 5 = 0.5), initially decreases until it 
reaches a minimum at two points symmetric with respect to x = 1/2, after which it starts 
increasing and keeps on doing so along the mixing-demixing line as the density is increased. 
This behavior of the critical temperature is related to the change in the order parameter of 
the transition. Close to x — or x — 1, where phase separation is essentially of liquid- vapor 
type, increasing the amount of the dilute component increases the weight of the interaction 
between unlike species in the internal energy. Since here we have 5 < 1, this leads to a 
decrease of the overall attractive contribution to the internal energy, resulting in a lower 
critical temperature of the liquid-vapor transition. On the other hand, the same argument 
implies that for a transition which is mainly mixing-demixing in character, approaching 
equimolar concentration increases the gain in internal energy entailed by the demixing, so 
that the critical temperature increases as one moves towards x = 1/2. Once the mixing- 
demixing critical line has been reached, an increase of the density at constant concentration 
similarly favors the energetic contribution to the free energy and leads to an increase of the 
critical temperature. 

As 5 is increased, the mixing-demixing line, as noted above, moves to higher density, 
and so does that portion of the line originating from the pure species where the transition 
is predominantly of demixing type. At about 5 = 0.46, a new feature appears in the critical 
lines, namely a crescent-shaped line at low density and concentration spanning an interval 
centered at x = 1/2, where the transition is essentially liquid- vapor. This is shown in Fig. |2| 
for 5 = 0.65. The occurrence of such a critical line can be euristically explained as follows: 
if 5 were equal to 1, the mixture would reduce to a one-component hard-sphere fluid with 
attractive tail interaction. If 5 is not too small, a nearly equimolar mixture may still behave 
like a sort of "effective" one-component fluid displaying a liquid-vapor transition. In order 
for this to happen, however, the density has to be low enough, so that the smaller gain 
in internal energy resulting from choosing a liquid-vapor phase separation at intermediate 
concentration instead of a mixing-demixing one can be compensated by a larger entropy. 
We must point out that, on the basis of the investigation performed in Ref. |lj and the 
discussion made here below Eqs. ([TT]) (p~5|) , we do not expect this critical line to appear in 
the equilibrium phase diagram right above 5 = 0.46. In fact, the interval 0.46 < 5 < 0.605 
corresponds to the hidden-binodal regime of Ref. where the critical point at x = 1/2 of 
the crescent-shaped line is metastable. Of course, knowledge of the behavior of the point 
at equimolar concentration alone is not sufficient to deliberate about the fate of the whole 
critical line. In principle, some portions of it might become stable for 5 different from the 
value 5 = 0.605 reported in JTJ. However, this would imply a change in the topology of the 
critical lines with respect to that shown in Fig. |2] which we did not observe when fluctuations 
are taken into account (see Sec. |TV| ), and we regard such an occurrence as rather unlikely. In 
summary, it is palusible that the value 5 = 0.605 obtained in Ref. [I] sets the threshold for 
the appearance in the mean-field equilibrium phase diagram of the whole crescent-shaped 
critical line besides the point at x = 1/2. The critical temperature along the crescent- 
shaped line changes very little, typically by few percent, and it presents a shallow minimum 
at x = 1/2. Along the other critical lines the qualitative behavior of the temperature is the 
same as in Fig. [I]. For relatively low S, the minima of the critical temperature located along 
the critical line arising from the pure species are higher than the temperature at the ends 
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of the crescent- shaped line, while the converse is true for S larger than about 0.64. 

As S is increased, the crescent grows towards larger and smaller concentrations, until for 
5 = 0.65338 the critical lines meet each other, resulting in the topology shown in Fig. |3|. 
When 5 grows above this value, the former critical line originating from the pure components 
splits: the portions at low and high concentration, where the transition is mainly liquid- 
vapor, join the crescent-shaped line so as to form a liquid-vapor critical line ranging from 
x = to x = 1, while the part at intermediate concentration, where demixing prevails, 
remains connected to the demixing line at x = 1/2 and detatches from the liquid- vapor 
line, giving a fork-shaped critical locus. The situation just described is illustrated in Fig. |] 
for 5 = 0.7. The liquid-vapor critical line has just one temperature minimum at x — 1/2, 
in agreement with the above observation that for a liquid-vapor transition, increasing the 
concentration of the dilute component leads to a decrease of the critical temperature. The 
two temperature minima symmetric with with respect to x = 1/2 are now located along the 
fork-shaped line. The relative temperature change along this line is however quite small, as 
already observed for the crescent-shaped line of Fig. |^. For 5 just above the value 0.65338 
that marks the boundary between the topology of Fig. |2| and that of Fig. |], the critical 
temperature at the tips of the fork-shaped line is higher than the minimum at x = 1/2 on 
the liquid-vapor line, while the converse is true at higher 5, including the value 5 = 0.7 to 
which Fig. ^ refers. A similar behavior is found when comparing the minimum on the liquid- 
vapor line with the temperature at the intersection of the fork with the mixing-demixing 
line. 

If S is further increased, the liquid- vapor line becomes more and more similar to a straight 
segment, as is to be expected since 5=1 corresponds to a one-component fluid, whose 
critical density is obviously independent of the concentration. At the same time, the fork- 
shaped line shrinks and moves to higher density, together with the mixing-demixing line. 
Strictly speaking, the fork disappears from the critical lines only in the one-component limit 
S — > 1. However, when 6 gets larger than a value 5o between 0.75 and 0.76, this locus is 
certainly metastable, as the pressure along it is everywhere negative. We observe that this 
is a sufficient condition for metastability, but not a necessary one. In fact, according to 
Ref. 0, the portion of the fork near equimolar concentration has already disappeared from 
the equilibrium phase diagram for 5 > 0.708. We are then left with two disconnected critical 
lines: the mixing-demixing one, that terminates at an endpoint, and the liquid- vapor one 
at lower density, as shown in Fig. [| for 5 = 0.8. For S — > 1, this topology evolves into 
that expected for the one-component fluid, as explained above: the mixing-demixing line 
eventually disappears, and the liquid-vapor line becomes a segment at constant density. 



IV. HRT PHASE DIAGRAM 

As we said in the Introduction, HRT calculations were performed for an interaction that 
consists of a hard-sphere repulsive core and an attractive Yukawa tail, whose inverse-range 
parameter z was set to z = 1.8 for both like and unlike species. This value of z is appropriate 
for representing the LJ potential |3T| and it has been widely adopted in the literature. We 



illustrate our results by considering, for different values of the parameter <5, several isothermal 
sections of the phase diagram, each of which corresponds to a single HRT run. In our opinion, 
this gives a clearer picture than the one that would be obtained by mapping the phase 
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diagram at given 5 on a single three-dimensional plot. In the following, the temperature 
will be identified with the corresponding reduced quantity T* = k B T/e, and the asterisk 
will be omitted. Fig. |6| shows the phase diagram on the density-concentration plane at 
four different temperatures for 5 = 0.65. Because of the symmetry of the model considered 
here, the phase diagram is obviously symmetric with respect to the equal-concentration 
axis x = 1/2. We already pointed out that in the HRT the conditions of thermodynamic 
equilibrium that define the coexistence region are implemented by the theory itself. In 
fact, inside the domains shown in the figures the hessian determinant of the Helmholtz 
free energy is identically vanishing. We stress that this is essentially different from what is 
found in mean field-like approaches, where phase separation is marked by the appearance 
of spinodal surfaces in the T-p-x space, which give spinodal lines upon intersecting with a 
plane at constant T like those of the figure. Inside the regions bounded by these lines, the 
hessian of the free energy attains unphysically negative values. In the HRT, on the other 
hand, the hessian does not become negative, but it vanishes identically in a region of finite 
|32jl . If we consider a curve at given temperature and pressure in the p-x plane, 
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it is readily seen that along the portion of this curve that lies inside the domain where the 
hessian vanishes, the chemical potential of both species are identically constant. Therefore, 
such a domain is indeed the coexistence region of the mixture. This clearly appears from 
Fig. |7], which shows that the domains of Fig. collapse into lines when they are plotted 
in the P-Ap plane, where P is the pressure, and Ap — pi — p 2 is the difference between 
the chemical potentials of the two species. Because of the conditions of thermodynamic 
equilibrium, the coexistence regions in the T-P-Ap space appear as "sheets" bounded by 
critical lines. Intersecting with a plane at constant T then yields lines like those of Fig. [7], 
terminating at critical points. Each point of these lines corresponds to an isobar of the 
domains of Fig. |6], i.e., to a tie-line. At a critical point, the tie-line reduces to a single point. 
We recall that in the p-x plane a critical point is not in general an extremal point of the 



phase boundary pD| , either in p or in x, and it cannot be detected by just considering the 
shape of the p-x coexistence boundaries. Critical points have been marked by dots in Fig. |6| 
and in the following figures that show the phase diagram in the p-x plane for different values 
of 5. 

The first panels of Figs. || [7| show the phase diagram at a temperature somewhat lower 
than the critical temperature T c ° ~ 1.2 of the pure components. Three distinct coexistence 
domains are present. The two at low density originate from the coexistence regions of the 
pure components, and they both terminate at a critical point. The high-density region at 
Ap = 0, which is present also at temperatures above T c °, involves, at each fixed density, 
coexistence between two fluids at concentration x and 1 — x respectively, and terminates 
at a mixing- demixing critical point at x — 1/2. Below a certain temperature T t ~ 1.06, 
however, the demixing region bifurcates into two branches (see the second panel of Fig. 0), 
each of which ends at a critical point with Ap ^ 0, x ^ 1/2. The former mixing-demixing 
critical point at x = 1/2 has now become a first-order coexistence boundary, at which a 
mixed fluid at equal species concentration coexists with a demixed fluid at higher density 
similar to that found above T t , consisting of two phases with concentrations symmetric 
with respect to x = 1/2. At the temperature T t at which the bifurcation develops, the 
critical line at x = 1/2 and the two critical lines generated by the symmetric branches for 
T < T t meet at a tricritical point. These critical lines can be visualized as the boundaries 
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of the coexistence "sheets" in the T-P-Ap space whose projections on the P-Ap plane at 
constant T are shown in Fig. [7|. At the tricritical point, the mixed fluid at x = 1/2 and 
the two phases that constitute the demixed fluid become critical simultaneously. Again, 
in the p-x plane the tricritical point does not display any special feature that makes it 
immediately detectable. Inspection of the p-x phase boundary alone does not allow one to 
tell whether the high-density coexistence region displays a single critical point at x = 1/2, 
two symmetric critical points, or a tricritical point. By further lowering the temperature, 
the high- and low-density coexistence regions expand and get closer, until at a temperature 
Td ~ 1.04 each of the low-density regions meets the high-density one at a double critical 
point. The critical lines in the p-x plane have the same topology as in Fig. [I]. In particular, 
the tricritical point corresponds to the intersection of the mixing-demixing critical line at 
x = 1/2 with the line that spans the concentration axis from x = to x = 1. The latter 
actually results from projecting on the p-x plane the two critical lines into which the mixing- 
demixing critical line bifurcates below T t and those originating from the critical points of 
the pure species. The two symmetric double critical points at T = where these critical 
lines meet in couples correspond to the temperature minima at x ^ 1/2 located on the 
critical line of Fig. [|. We observe that, unlike tricritical points, double critical points do not 
entail the intersection of topologically distinct critical lines. In fact, tricritical points are 
found in the model considered here because of its special symmetry, but they do not occur 
in real binary mixtures, while double critical points are frequently found in real mixtures, 
including mixtures of noble gases such as neon- krypton and neon-xenon [|3~4}| . Below Td, 
the coexistence domain consists of one connected region, without any critical point, as shown 
in the third and fourth panels of Fig. 0. A section of the phase diagram at x = 1/2 in the p-T 
plane is plotted in Fig. |||, showing the same behavior as in Fig. 2(d) of Ref. [|TJ. The HRT 
and the mean- field critical lines are compared in Fig. ^ for 5 = 0.5. This value of S is small 
enough to give the same topology of the critical lines in HRT and mean-field theory, save for 
an extremely short line at low density and nearly equimolar concentration obtained in mean 
field, which belongs to the metastable regime (see the discussion in the previous Section) 
and has not been reported here. The figure shows the projections of the critical lines both in 
the p-x and in the x-T planes; in the latter case, the mixing-demixing line at x — 1/2 has not 
been shown. For the pure species, it is known from the comparison of the mean-field results 
with accurate simulation data for the critical constants of the Yukawa fluid with the same 



inverse range z = 1.8 considered here [pq| , that mean-field theory underestimates the critical 
density by about 20% and overestimates the critical temperature by about 10%. Basically 
the same differences are found by comparing the mean-field and the HRT results, as the HRT 
provides a very good determination of the critical point of LJ-like fluids P^,55,57||. Fig. |^ 
shows that similar discrepancies between mean field and HRT hold also for the critical loci 
of the binary system. We observe that the direction of the order parameter is little affected 
by the inclusion of fluctuations, at least for the present case where these do not change the 
topology of the critical lines. 

According to our HRT calculations, the value 5 = 0.65 is actually nearly coincident with 
the upper limit of 5 for the topology that we have just described. The scenario at slightly 



larger 5 is illustrated in Figs. [TO, O, which show the phase diagram in the p-x and P-Ap 



planes for 5 = 0.665. At temperatures T above 1.023, the phase diagram evolves as before: 
on lowering T, the high-density coexistence region bifurcates at a tricritical point, and each 
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of the resulting branches merges with the low-density coexistence regions that originate 
from the pure species, so that just below T = 1.03 there are no critical points left in the 
phase diagram. However, at T = 1.023 two critical points reappear at low density. In the 
P-Ap plane, this is marked by the appearance of two "twigs" that stick out of the low- 
density coexistence region, each terminating at a critical point, as shown in the third panel 
of Fig. [□]. By slightly lowering T, these twigs grow longer, and for T = T d = 1.022 they 
meet at a critical double point located at Ap = 0. At the same time, the low-density lobes 
of the coexistence region in the p-x plane become very elongated, and coalesce at x — 1/2. 
The topology of the critical lines is the same as in Fig. where the crescent-shaped line at 
low density is that described by the new family of critical points. We remark that, according 
to mean-field theory, this line is extremely shallow with respect to the temperature. This 
is confirmed by the HRT results just reported, which give a relative variation of the critical 
temperature of about 0.1%. Below T d (see the fourth panel of Fig. |10]) , the phase diagram 
presents an "island" of homogeneous, mixed fluid surrounded by a "sea" of phase-separated 
fluid. A section of the coexistence region at x = 1/2 shows that in this regime we have two 
kinds of phase equilibria at equal species concentration: on the one hand, below the tricritical 
temperature T t a demixed fluid at high density coexists with a mixed one at intermediate 
density as before. On the other hand, this mixed fluid coexists, at slightly lower density, 
with another low-density mixed fluid. The latter phase equilibrium terminates for T = T d 
at a vapor-liquid critical point. Therefore, the phase diagram at x — 1/2 exhibits both a 
critical and a tricritical point. On further lowering the temperature below T d , the island of 
mixed fluid shown in Fig. [X0| becomes smaller and smaller, until it is eventually swallowed 
by the coexistence region. Correspondingly, in the p-T plane the two phase boundaries at 
x = 1/2 meet at a triple point, where the two mixed fluids at different densities coexist 



with the demixed fluid P5[ . Below the temperature of the triple point, the demixed fluid 
at high density coexists with the mixed one at low density. The phase diagram at equal 
species concentration is similar to Fig. 2(c) of Ref. |TJ and to that shown below in Fig |TJ| 
for the case S = 0.67. From the picture given above it is clear that the low-density lobes of 
the coexistence region meet at a temperature lower than that at which each lobe meets the 
high-density region. This in turn is lower than the temperature at which the tricritical point 
appears. Therefore, for the topology just discussed the liquid-vapor critical temperature T d 
is always lower than the tricritical temperature T t . 

A different situation arises as 5 grows above 0.67. Figs. 13, [14] depict the evolution of 
the phase diagram for 5 = 0.68. As before, a tricritical point is present at x = 1/2, as 
shown by the bifurcation of the high-density coexistence region. On the other hand, on 
lowering the temperature the low-density branches of the coexistence region do not merge 
with the high-density one as before, but instead meet each other at a critical double point 
between T = 1.03 and T = 1.025. In this regime the isothermal sections of the phase 
diagram in the p-x plane consist of two disconnected domains, with that at lower density 
spanning the whole density-concentration axis from x — to x — 1. At about T = 1.025, 
two twigs sprout out of the low-density coexistence domain in the P-Ap plane. In such 
a situation this domain presents two symmetric critical points, corresponding to the ends 
of the twigs. On further lowering T, the twigs grow longer, and for T between 1.015 and 
1.01 each of them meets a branch of the high-density coexistence region at a critical double 
point. The low- and high-density coexistence domains then join at two points symmetric 
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with respect to x = 1/2, leaving between them an island of one-phase fluid (see Fig. |13| ), 
that eventually disappears at low temperature. This scenario is somehow the converse of 
that previously described for 5 = 0.665, where the high-density coexistence domain merges 
with the low-density ones originating from the pure components, and the one-phase island 
inside the coexistence region is formed because of the twigs meeting each other at x — 1/2. 
The topology of the critical lines for the case just discussed is that of Fig. |j: this time the 
critical line that goes from x = to x = 1 is not connected to the mixing- demixing critical 
line, and it has a temperature minimum at x — 1/2, corresponding to the coalescence of the 
low-density coexistence regions into one connected domain. The mixing-demixing critical 
line intersects at the tricritical point two symmetric critical lines as before, which however 
do not meet those originating from the pure species. Instead, they merge with the lines 
described by the family of critical points that sprout out of the low-density coexistence 
region, resulting in the fork-shaped line of Fig. f|. This presents two symmetric temperature 
minima at the value of T at which the critical lines meet each other. As already observed for 
the crescent-shaped line of Fig. 0, the temperature along the fork is actually almost constant, 
both according to mean-field theory and HRT, so that the minima are extremely shallow. 
The topology of the phase diagram at equimolar concentration is similar to the previous one 
illustrated in Fig. [T2| for 5 = 0.67. In particular, in a certain temperature interval including 
that where the two-phase region encloses a one-phase domain, there are two different phase 
equilibria at x — 1/2. For the coupling parameter 5 = 0.68 considered here, the liquid- vapor 
critical temperature Td, which is the temperature minimum along the liquid-vapor critical 
line, is smaller than the tricritical temperature T t , as for 5 = 0.665. However, in the present 
topology where, on lowering T, the low-density coexistence domains merge with each other 
before they merge with the high-density one, nothing prevents the liquid- vapor critical point 
at x = 1/2 from occurring at a higher temperature than the tricritical point. In the P-Ap 
plane, this corresponds to the high-density coexistence region bifurcating at a temperature 
lower than that at which the low-density branches meet. Such a scenario in fact comes along 
starting from about 5 = 0.7. The phase diagram at x — 1/2 for this value of 5 is shown in 
Fig. 15, and is the same at that of Fig. 2(b) of Ref. Q. 

The phase behavior of the system for 5 intermediate between the values 0.665 and 0.68 
discussed above is worth being considered in more detail. Figs. [If], [17] show the phase diagram 
for S = 0.67. On lowering the temperature, the low-density coexistence domains coalesce at 
about T = 1.0235, and at a slightly lower temperature around T = 1.023 the resulting region 
coalesce with the high-density coexistence domain. The phase diagram in the p-x plane is 
then similar to that already shown for 5 = 0.68, except that here the low- and high-density 
domains meet at nearly the same temperature. However, if one considers temperatures just 
above T = 1.0235 it appears that, unlike what found for the 6 studied above, the phase 
diagram displays six critical points, two for each of the disconnected domains that make up 
the coexisting region. This is clearly shown in the P-Ap plane of Fig. [T7| for T = 1.024. 
In the present case it is not obvious to tell whether the topology of the critical lines in 
the p-x plane is that of Fig. ^| or of Fig. ^, because it is difficult to ascertain which is the 
part of the coexistence region in the P-Ap plane that is sprouting from the other, and 
consequently whether the high-density branches of the coexistence domain are bound to 
meet those originating from the pure species as for S = 0.655, or the twigs as for S = 0.68. 
We are therefore very close to the boundary between the two topologies of the critical loci, 
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that corresponds to the scenario of Fig. |3|. This is marked by the presence of two tricritical 
points at non-equimolar concentration besides that at x = 1/2. In such a situation the 
growth of a secondary structure (the twigs) from an already existing one (the main branch) 
is replaced by a bifurcation where both branches stem out of the tricritical point at the 
same temperature. For continuity reasons, we expect that in a very narrow range of 5 values 
that includes 0.67 and is contained in the interval 0.665 < 5 < 0.68, the configuration with 
six critical points shown in Fig. [T7] can be found for the topology of the critical lines of 
both Fig. |2| and Fig. |j. If one considers the critical lines of Fig. 0, this occurs when the 
critical temperature at the ends of the low-density crescent-shaped line is higher than the 
minima located along the line originating from the pure species. In the case of Fig. f|, the 
requirement is that the local temperature minimum at x = 1/2 along the liquid- vapor line 
must be lower not only than the tricritical temperature, but also than that of the critical 
points located at the tips of the fork-shaped line. It may also be worth pointing out that, 
strictly speaking, the watershed between the density-concentration phase diagrams of Fig. |H] 
and Fig. [13], that represent the two different ways by which a domain of one-phase fluid can 
be enclosed into the coexistence region, has not to coincide with the boundary between the 
critical lines just discussed, although the two are expected to occur for very similar values of 
5. In fact, according to mean-field theory, when the critical lines have the topology of Fig. |3|, 



the coexistence region in the p-x plane still looks like that of Fig. |10| This implies that 
there is an extremly narrow range of 5 where the topology of the critical lines is as in Fig. |], 
and the topology of the phase diagram in the p-x plane is as in Fig. [U]. In the P-A/i plane 
such a regime is marked by the twigs joining the bifurcation of the high-density coexistence 
region at a temperature higher than that at which the low-density branches meet each other. 
We have not checked whether this scenario comes along also in the HRT, or instead it is 
replaced by the converse one, where the critical lines have the topology of Fig. |2], and the 
phase diagram in the p-x plane has the topology of Fig. [T~3| . 

Further information on the nature of the phase equilibria is obtained by considering 
the tie-lines, i.e., the lines in the p-x plane that connect the phases at coexistence at a 
certain temperature and pressure. A few tie-lines of the mixture with 5 = 0.67 are shown in 
Fig. [18| for T = 1.024 (left panel) and T = 1.023 (right panel), corresponding respectively 
to the second and fourth panel of Fig. [16]. We have chosen to report the tie-lines for the 
case 5 = 0.67 because of its particularly rich phase diagram, several features of which are 
separately found also for different values of 8. In both panels of Fig. [Tj| the high-density 
portion of the coexistence region is characterized by the presence of a demixed fluid which, 
as said above, consists of two phases at the same density and concentrations symmetric with 
respect to the equal-concentration axis. As a consequence, the average density coincides 
with that of the coexisting phases irrespective of their relative amount in the demixed fluid, 
resulting in strictly vertical tie-lines. Since both panels refer to temperatures below the 
tricritical temperature T t , the phase equilibrium just described does not extend down to 
the left boundary of the high-density coexistence domain. At a certain density p-£>, the two 
coexisting fluids at symmetric concentrations x, 1 — x coexist in turn with a fluid at a lower 
density pm and equimolar concentration. The densities pm, Pd are the boundaries of the 
first-order coexistence domain that separates the mixed-fluid and demixed-fluid regions of 
Fig. [12] for T <T t . In the p-x plane, this three-phase equilibrium takes place in a triangle- 
like domain bounded by three tie-lines meeting in couples at the points (pn,x), (pd, 1 — 
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(Pm, 1/2). These lie very close to (but do not exactly coincide with) the two oblique lines and 
the leftmost vertical line shown in the high-density region of Fig. |18[ The two symmetric 
portions of the high-density region of the left panel that lie above and below the three-phase 
domain correspond in the P-Ap plane to the branches of the bifurcation that stems from 
the line at Ap = (see Fig. |I7|). Each of them presents coexistence between two phases 
which differ both in density and concentration, terminating at a critical point where the 
tie-line reduces to a single point. In the right panel the critical points are absent, as the 
high- and low-density coexistence regions have merged into one connected domain. The 
left panel also shows the tie-lines shrinking at the four critical points located on the two 
symmetric lobes of the low-density coexistence region. An interesting feature that appears 
from the arrangement of the tie-lines in this region is the presence of two more domains of 
three-phase coexistence, such that the coexisting phases have concentrations that lie on the 
same side of the x — 1/2 axis. At temperatures at which the coexistence region is connected 
and encloses an island of homogeneous fluid as in the right panel, this sort of three-phase 
equilibrium could be expected just on the basis of the fact that different tie-lines cannot 
intersect each other. However, the presence in the low-density region of the two critical 
points associated with the twigs in the P-A/i plane implies that three-phase coexistence is 
observed in this region even at temperatures at which it has not yet merged with the high- 
density coexistence domain. We also observe that in the neighborhood of the x — 1/2 axis 
where the two lobes of the low-density region coalesce, the tie-lines are nearly horizontal, 
meaning that the concentration of the coexisting phases are very similar. At exactly x = 1/2 
azeotropy occurs, i.e., one finds a purely liquid- vapor transition between two phases at equal 
concentration, as already observed above in connection with the phase diagram at x = 1/2 
in the p-T plane. The critical point topping the liquid-vapor coexistence curve at equimolar 
concentration in Figs. O, [TBI is therefore an azeotropic critical point. 

As S increases above 0.7, the curvature of both the low- and high-density coexistence 
boundaries in the p-x plane rapidly decreases, and the pressure at which the bifurcation of the 
high-density coexistence region takes place gets closer and closer to that of the coexisting 
vapor and liquid phases on the low-density region at x — 1/2. As a consequence, the 
one-phase island bounded by the coexistence region becomes smaller and involves a fluid 
at nearly equimolar concentration. One can then ask whether this feature of the phase 
diagram will disappear at a certain Sq < 1, or instead only in the limit 8 — > 1. The former 
case corresponds to a situation where for 5 > 5o the high-density coexistence region does 
not undergo any bifurcation in the P-Ap plane, but it always presents a mixing-demixing 
critical point at x — 1/2, Ap = 0. At low enough temperature, this critical point meets 
the low-density coexistence region at a first-order phase boundary. The topology of the 
critical lines in the p-x plane is that of Fig. |5|: the mixing-demixing line terminates at 
a critical enpoint where the low- and high-density coexistence regions meet. The phase 
diagram at x = 1/2 corresponding to this scenario is that of Fig. 2(a) of Ref. JTJ] . As 
before, the equimolar mixture displays both a liquid-vapor transition that occurs below 
a certain critical temperature and involves two mixed fluids of different densities, and a 
mixing-demixing transition at high density. In this case, however, the transition between 
the mixed liquid and the demixed fluid is always second-order down to the temperature of 
the endpoint, below which the demixed fluid coexists with a mixed vapor of much lower 
density. The situation where the contact between the low- and the high-density coexistence 
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region always occurs at two distinct points symmetric with respect to x = 1/2 until S gets 
equal to 1 is instead consistent with the phase diagram remaining qualitatively similar to 
that found for 5 = 0.7 with critical lines as in Fig. f|, albeit the fork-shaped line will become 
vanishingly small as 5 approaches 1 and the mixing-demixing line will be pushed to higher 
density. The phase diagram of the equimolar mixture in the p-T plane will qualitatively 
look as that of Fig. [HJ As a consequence, as long as one has S < 1 there will be a small 
temperature interval between the triple point and the tricritical point, where the transition 
between the mixed and the demixed liquid is first order. This is the scenario advocated in 
Rcf. PU[ on the basis of modified hypernetted chain (MHNC) calculations on a symmetric 



LJ mixture, while the simulations performed in Ref. [JXJ] on a square-well system supported 
the existence of a critical endpoint at high enough 5. Some calculations that we previously 
performed on a hard-sphere LJ mixture using a 100 x 100 density-concentration grid also 
suggested that for 5 = 0.8 the phase diagram of the system presents a critical endpoint. 
However, the present results for a hard-core Yukawa potential shown in Figs. [T^, PU| indicate 
that for 5 = 0.8, the mixture most likely does not have a critical endpoint: in fact, in a certain 
narrow temperature range the coexistence region in the p-x plane still encloses a domain of 
one-phase fluid with concentration varying in a small interval around x — 1/2. Since the 
inverse range z = 1.8 of the Yukawa potential considered here gives a fair representation of 



the LJ interaction [31], the difference between this behavior and that reported in Ref. [19 
most probably does not depend on some intrinsic difference between the two interactions, 
but just on the fact that the higher resolution entailed by the 150 x 150 grid used in this 
work allows one to study the contact between the low- and high-density coexistence regions 
with better accuracy than in Ref. |19| , thereby uncovering the presence of a small one-phase 
domain. We should observe that the case 6 = 0.8 actually lies at the limit of resolution even 
for the larger grid employed here: the homogeneous domain inside the coexistence region 
in the p-x plane extends only one gridpoint in the p-direction, while the phase diagram 
in the P-Ap plane does not present any of the features associated with the occurrence of 
such a domain on the scale of the figure, so that it cannot be clearly distinguished from 
the kind of diagram one would expect in the presence of a critical endpoint. In particular, 
there is no sign of the bifurcation of the high-density coexistence region in the P-Ap plane 
that we expect if the low- and high-density regions have to meet at two points of non- 



equimolar concentration according to the scenario of Fig. [L4l This most probably depends 



on the fact that, for the relatively high value of 5 considered here, the densities of the mixed 
and demixed liquids at coexistence are so close, that their difference is comparable to the 
mesh size. As a consequence, the differences in the pressures and chemical potentials of 
neighboring but non-coexisting points will also become comparable to the numerical errors. 
A little numerical noise in the P-Ap plane due to the finite mesh used in the calculation 
is in fact always present, as shown by close inspection of the relevant figures, but only at 
high 5 does this become a hindrance for a clear description of the phase behavior. We did 
not pursue any investigation for 5 > 0.8, because a fortiori in such a regime the resolution 
allowed by the grid used here would not enable us to discriminate between the occurrence 
of a critical endpoint, or an extremely small region of homogeneous fluid that undergoes a 
first-order demixing transition on slightly increasing the density. Therefore, we are not in 
a position to say whether a critical endpoint at x = 1/2 will eventually appear near 6=1, 
or instead a very weak first-order transition will survive up to 5 = 1. Another issue that 
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is difficult to elucidate within the resolution of the present calculation concerns the fate of 
the twigs discussed before in connection with the phase diagram topology of Figs. 13 , 14. 
As noted above, even at temperatures at which the low-density coexistence regions that 
originate from the pure species have coalesced into one connected domain, the latter can 
still show two symmetric critical points. When this is the case, the low-density coexistence 
region in the P-Ap plane presents two small twigs, each of them ending at a critical point. 
These twigs could either disappear at high enough 5, or persist for 5 arbitrarily close to 
1. If a critical endpoint at x — 1/2 never appears and the topology of the critical lines in 
the p-x plane is that of Fig. |] until S gets equal to unity, the latter case corresponds to a 
situation where the low- and high-density coexistence regions will always meet at two double 
critical points. As discussed above, mean-field theory indeed gives two double critical points 
located along the fork-shaped line, but for 5 close to 1 they are certainly metastable since 
they correspond to negative pressure. If instead the twigs disappear before 5 reaches 1, there 
will be a value of S above which the critical points of the high-density coexistence region meet 
the low-density region at two points of first-order coexistence. As a consequence, the contact 
takes place at two critical endpoints with concentrations symmetric with respect to x = 1/2 
rather than at two double critical points. This is the situation considered in Ref. P5|| . In 
such a case the critical lines in the p-x plane have still the same topology as in Fig. [|, except 
that the temperature minima along the fork-shaped line have now to coincide with the tips of 
the fork, which correspond to the two endpoints. It can also be worthwhile observing that in 
principle the twigs are not incompatible with a critical endpoint at equimolar concentration. 
In fact, this could even give rise to a scenario where the presence of the endpoint is consistent 
with the high- and low-density coexistence regions meeting at non-equimolar concentration. 
Specifically, if the critical points located at the ends of the twigs meet the high-density 
coexistence region at two symmetric critical endpoints with x ^ 1/2, nothing prevents the 
point at equimolar concentration of this region from remaining critical, until it also meets 
the low-density region at a critical endpoint. We remark that such a possibility hinges on the 
presence of the twigs and the related critical points on the low-density coexistence region. 
In the converse situation where the critical endpoints at x ^ 1/2 result from the contact 
between two points of first-order coexistence on the low-density domain and two critical 
points on the high-density one, the requirement that the point at x = 1/2 be critical is 



untenable, as pointed out in Ref. p9fl . The scenario just depicted is somewhat suggestive, as 
it could account for both the presence of a critical endpoint at x = 1/2 above a certain value 
of S reported in simulation studies M, and the failure to observe the low- and high-density 



coexistence regions coalescing at x = 1 /2 found in the present investigation as well as in [|39| 



However, we must point out that at the present stage we do not have any solid evidence 
that this possibility does actually occur, so that we must regard it as a purely speculative 
conjecture. 



V. CONCLUSIONS 

We have used the HRT to perform an investigation of the phase diagram of symmetric 
binary mixtures as a function of the unlike-to-like interaction ratio 5. The microscopic in- 
teraction adopted consisted of a hard-core repulsion plus an attractive Yukawa tail potential 
with inverse range z = 1.8. Such a potential has been used many times in liquid-state 
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theory to describe a simple LJ-like fluid. For each value of S considered, results for the 
coexistence regions were obtained on the whole density-concentration plane at several tem- 
peratures. The resulting phase portrait was related both to the different topologies of the 
mean-field critical lines that come along as 6 is varied in the interval < 5 < 1, and to the 
behavior predicted by mean-field theory and simulation results [0] for a certain special class 
of the systems considered here, namely equimolar mixtures at molar fraction x — 1/2. A 
feature of the HRT that is particularly useful is that the domains of coexisting phases are 
straightforwardly obtained by the theory as the loci where the conditions of thermodynamic 
equilibrium between different phases are satisfied, without any need of enforcing them a 
posteriori. Fulfillment of these conditions is shown by the collapse of the isothermal sections 
of the coexistence regions on lines of the P-A/i plane, P being the pressure, and A/z the 
difference between the chemical potentials of the two species. This property makes also easy 
to identify the critical points exhibited by the phase diagram at a certain temperature. 

According to the results previously obtained by mean-field theory and simulations, the 
phase diagram of symmetric mixtures at equimolar concentration is characterized by three 
different regimes, depending on the value of 5: at low 5 (5 < 0.605 according to mean- 
field theory |]]) a mixing-demixing critical line joins a mixing-demixing coexistence curve 
at a critical point, thereby generating a tricritical point. For 5 closer to 1 (5 > 0.708 in 
mean field [|l|), the mixing-demixing critical line joins a liquid- vapor coexistence curve at a 
critical endpoint. Finally, in a narrow interval in 5 intermediate between the above regimes, 
the mixture shows both a liquid-vapor coexistence curve and a mixing-demixing one at 
higher densities, topped respectively by a critical and a tricritical point. Our investigation 
clearly shows the tricritical point regime as well as the intermediate one. The latter is 
predicted to occur starting from about 5 > 0.65, in agreement with the simulation results 
for a square- well mixture [Q], but it lingers on for larger values of 5 than those given by 
simulation: according to simulation, the endpoint regime is reached for 5 > 0.68, while 
according to HRT 5 = 0.7 is still in the "transition" or intermediate regime, as shown in 



Fig. [15]. As a matter of fact, we did not find any clear evidence of the endpoint regime, since 
our calculations indicate that the mixture is likely to be in the intermediate regime for 5 as 
high as 0.8. However, on increasing 5, the transition region quickly moves to high density, 
and the finite resolution allowed by the 150 x 150 grid used here becomes insufficient to 
fully uncover the topology of the phase diagram in this parameter range. For this reason 
we did not investigate the behavior of the model for 5 > 0.8. The results obtained in this 
work are in qualitative agreement with a study based on the MHNC integral equation for 



a LJ mixture |j39| , according to which no critical endpoint at equimolar concentration is 



present up to at least 5 = 0.81. Before drawing any definite conclusion about the resilience 
of the intermediate topology in the theoretical HRT and MHNC results compared with the 
simulations and the mean-field approximation, the role of the specific interaction adopted 
should be elucidated. In fact, the HRT and MHNC calculations were performed on a HCY 
and on a LJ potential respectively, while in the simulations a hard-sphere plus square-well 
was used. As observed in Sec. Q, the HCY potential was also employed in other investigations 
based on the MSA J0J, the ORPA g2f, and the SCOZA |3| theories, all of which yield for 



the phase diagram the same qualitative picture found in mean-field theory. However, we do 
not see any general reason why the independency of the interaction profile which is intrinsic 
to the mean-field phase behavior should always hold also for more sophisticated approaches, 
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specifically for the HRT. Moreover, one can also imagine a scenario where the existence of a 
critical endpoint at equimolar concentration is not prevented by the behavior observed here 
for the phase diagram in the p-x plane, namely the fact that, even for relatively high 8, the 
low- and high-density coexistence regions meet first at two points symmetric with respect to 
x = 1/2 rather than at x — 1/2. However, the present limits in resolution of our numerical 
calculation do not allow us to say whether this conjecture is actually relevant for the system 
studied. 

An interesting issue that was considered here is how the phase portraits corresponding 
to the topologies mentioned above look like, if one moves off the plane of equimolar con- 
centration. In particular, the intermediate regime is marked by the presence, in a certain 
temperature interval, of a phase coexistence domain with a hole of homogeneous fluid inside 
it. This comes along according to two distinct scenarios that can occur as the temperature 
is lowered: for 0.65 < 8 < 0.67, the two low-density coexistence regions originating from the 
pure species meet the high-density coexistence region associated with the mixing-demixing 
transition, and subsequentely they meet each other, leaving a domain of mixed fluid enclosed 
inside the two-phase region. For 8 > 0.67, first the two low-density coexistence regions meet 
each other, and then the resulting connected domain meets the high-density coexistence re- 
gion at two points symmetric with respect to concentration x — 1/2, leaving some one-phase 
fluid in between. In the former case, the temperature T c of the liquid-vapor critical point at 
x = 1/2 is always lower than the tricritical temperature T t , while in the latter case T c can 
be either lower or higher than T t , depending on 8. The present calculation gives T t < T c 
starting from about 8 = 0.7. 

Because of the existence of both a liquid-vapor and a demixing transition, the critical 
loci show both a line that spans the concentration axis and connects the critical points 
of the pure components, and a mixing-demixing line at x — 1/2. The tricritical point 
topology corresponds to a situation where these lines are connected, and the character of 
the transition changes continuously from liquid-vapor at x = or x = 1 to mixing-demixing 
at x — 1/2. In the endpoint regime, which as said above was not observed in the present 
investigation, the critical lines are instead disconnected, and the transition along the line 
that connects the critical points of the pure species is essentially of liquid-vapor type. The 
intermediate regime between these topologies is characterized by the presence of a further 
critical line at concentrations ranging in a certain interval centered at x = 1/2, which 
was referred to above as either the "crescent" or the "fork" line. The number of critical 
points that are found for a certain isothermal section of the phase diagram depends on the 
relative location of the temperature extrema along the critical lines. The behavior found by 
the present HRT calculation agrees qualitatively with that given by mean-field theory. In 
particular, the system is predicted to have up to six critical points at a certain temperature. 
According to HRT, this occurs for a very narrow range of 8 values contained in the interval 
0.665 < 8 < 0.68. At the boundary between the "crescent" and the "fork" topologies of the 
critical lines, the system has two tricritical points at symmetric non-equimolar concentrations 
besides that at concentration x — 1/2. In HRT, this particular topology occurs for a value of 
8 close to 0.67, to be compared with the mean-field result 8 = 0.65338. In both cases, these 
values lie very near the boundary between the two different types of one-phase "holes" in the 
phase coexistence region described above. When 8 is such that mean-field theory and HRT 
predict the same qualitative topology of the phase diagram, the quantitative discrepancy 
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between the critical loci given by the two approaches is similar to that found for the pure 
species. 

This investigation shows that symmetric mixtures, despite their conceptual simplicity, 
exhibit a very rich phase behavior. Like real mixtures, these systems have both liquid- 
vapor and mixing-demixing transitions. However, while in real mixtures the liquid-vapor 
and the mixing-demixing regimes generally correspond to states that differ widely in density 
and pressure, in symmetric mixtures these transitions can instead be located in the same 
region of the thermodynamic space. As a consequence, they tend to compete with each 
other, so that even a small variation in the relative strength of the interactions expressed 
by the parameter 5 is sufficient to bring about significant qualitative changes in the phase 
diagram. The results presented here are also relevant for Ising ferrofluids in the presence of a 
magnetic field. More generally, they show that HRT is capable of providing a comprehensive 
description and resolving even subtle features of the phase behavior of the model. This 
ability could prove useful to systematically study beyond the mean-field level also the phase 
diagram of more realistic, non-symmetric model mixtures that depend one more than just 
one parameter. 
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FIGURES 



FIG. 1. Density-concentration projection of the mean-field critical lines of a symmetric mixture 
with interaction parameter 5 = eyij^-w = 0.4. The total density p is the sum of the densities p\, p2 
of the components, and the concentration x is defined as pi/ p. The open dots mark the locations 
of the local minima in the critical temperature. The arrows indicate the direction of the order 
parameter (see text). 

FIG. 2. Same as Fig. [I] for 5 = 0.65. 
FIG. 3. Same as Fig. |l| for 5 = 0.65338. 
FIG. 4. Same as Fig. [I] for 5 = 0.7. 

FIG. 5. Same as Fig. [l] for 5 = 0.8. 

FIG. 6. Isothermal sections of the coexistence region of a symmetric HCY mixture in the p-x 
plane according to the HRT. The inverse interaction range is equal to z = 1.8 and the interaction 
ratio is equal to 5 = 0.65. The dots mark the locations of the critical points. 

FIG. 7. Same as Fig. || in the pressure-chemical potential plane. Ap = p\ — P2 is the difference 
between the chemical potentials of the components. Note how the points of the coexistence regions 
shown in Fig. ^ collapse into lines as a consequence of the conditions of thermodynamic equlibrium. 

FIG. 8. Phase diagram in the p-T plane of the HCY mixture with 5 = 0.65 for the special case 
of equimolar concentration x = 1/2. Open and full dots denote respectively the A-line and the 
first-order phase boundary. 

FIG. 9. Critical lines in the p-x plane (upper panel) and x-T plane (lower panel) of the HCY 
mixture with 5 = 0.5. Dotted lines: mean-field theory. Full dots: HRT. In the upper panel, the solid 
line is a guide for the eye obtained by smoothly interpolating between the dots, and the arrows 
indicate the direction of the order parameter. The order parameter along the mixing-demixing 
critical line at x = 1/2 is not shown here for clarity and is parallel to the x-axis according to both 
mean field and HRT. In the lower panel, the mixing-demixing critical line is not shown. 

FIG. 10. Same as Fig. § for 5 = 0.665. 

FIG. 11. Same as Fig. for 5 = 0.665. 
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FIG. 12. Same as Fig. || for 5 = 0.67. The letters V, L denote respectively the vapor and mixed 
liquid phase. 

FIG. 13. Same as Fig. | for S = 0.68. 
FIG. 14. Same as Fig. |7| for 5 = 0.68. 



FIG. 15. Same as Fig. [TJ for 5 = 0.7. 
FIG. 16. Same as Fig. | for 5 = 0.67. 

FIG. 17. Same as Fig. [7| for 5 = 0.67. 

FIG. 18. Coexistence region in the p-x plane for 5 = 0.67 at T = 1.024 (left panel) and 
T = 1.023 (right panel), corresponding to the second and fourth panels of Fig. [l6|, showing a few 
tie-lines connecting phases at coexistence. The different shades of gray give a measure of pressure 
(black: low pressure; white: high pressure). 

FIG. 19. Same as Fig. ^ for 5 = 0.8. The critical points in the second panel are not shown, 
because for this value of 5 the resolution of our numerical calculation does not allow us to locate 
them. 

FIG. 20. Same as Fig. fj] for 5 = 0.8. In the first panel, the high-density coexistence region at 
A/j, = does not appear as its pressures are outside the scale of the figure. 
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